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We describe an implementation of Grover's fixed-point quantum search algorithm on a nuclear 
magnetic resonance (NMR) quantum computer, searching for either one or two matching items in 
an unsorted database of four items. In this new algorithm the target state (an equally weighted 
superposition of the matching states) is a fixed point of the recursive search operator, and so 
the algorithm always moves towards the desired state. The effects of systematic errors in the 
implementation are briefly explored. 
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I. INTRODUCTION 

Grover's quantum search 0, 0|_is one of the key algo- 
rithms in quantum computation [2j , allowing an unstruc- 
tured database to be searched more efficiently than can 
be achieved by any classical algorithm. It is most simply 
described in terms of a binary function /, with an n bit 
input register, permitting N = 2™ inputs, and a single 
output bit. This function has the value 1 for k desired 
{matching or satisfying) inputs and zero for the other 
N — k inputs, and can only be investigated by an oracle 
which returns the value for any desired input. Grover's 
original algorithm uses a quantum oracle which performs 
the transformation 



Is) — (_!)/<*) Is) 



(1) 



which is applied alternately with an amplitude amplifi- 
cation operator. Beginning from the state H n \0), where 
H n is the n qubit Hadamard and |j) indicates an n qubit 
quantum register containing the number j, the system is 
rotated towards an equally weighted superposition of the 
satisfying inputs. A measurement of the register after 
0(^/N/k) steps will return one of the satisfying inputs 
with high probability, while a classical search will take 
O(Nfk) queries on average. 

Grover's search is known to be optimal j|| when the 
number of matching inputs is known, but problems oc- 
cur when k is unknown. As the basic procedure is a 
rotation, once the desired state is reached further itera- 
tions will drive the system away from this state. Thus it 
is necessary either to estimate the value of k (for exam- 
ple, by approximate quantum counting 0,0) or to use 
an algorithm which is more robust to errors in the value 
oik. 

Recently Grovcr has described a new quantum algo- 
rithm 7], which overcomes this problem by driving the 
system asymptotically towards the desired state: the al- 



gorithm will always move towards the target and can- 
not overshoot. This might seem impossible, as unitarity 
means that any iterative algorithm cannot have a fixed 
point, but the new algorithm overcomes this by using 
a process which is recursive rather than iterative, and 
so the target state can act as a fixed point. A further 
consequence of the fixed-point behavior is that the new 
algorithm should be relatively robust to certain types of 
systematic error in its implementation. 



II. THEORY 

Grover's algorithm comes in many forms and here 
we describe just one of these, before relating it to the 
new algorithm 0. Consider the transformation 



URolftRfU\0) 



(2) 



where Rf is a phase oracle which applies a phase of <j> to 
all basis states satisfying the function / and i?o applies 
this phase shift to the initial state |0). If we take U = H n 
and 4> = 7T then this corresponds to the first iteration 
of the original Grover algorithm. Subsequent steps are 
obtained by applying the last four operations r times, 
thus applying successive rotations. Larger rotations can 
be defined using a recursive approach, by taking 



Vr+i = V r RoVjRfV r 



(3) 
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with Vq = U. For the original search algorithm IP = 
U and R' = R, and so each recursive operator simply 
corresponds to one of the iterative operators. 

Grover's new fixed point quantum search differs signif- 
icantly from this by choosing 4> — ir/3, so that R^ does 
not equal R. Thus the recursive operators are not sim- 
ply iterative, and have to be worked out separately for 
each value of r; for more details see the original paper 
Q. Here we consider the case of n — 2, with either k = 1 
or k = 2, and take U as the pseudo-Hadamard gate (a 
90° rotation), as these are the examples we implement 
experimentally. 
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TABLE I: Success probabilities, P r , and the total number 
of queries used, Q r , for the r th stage of Grover's fixed-point 
quantum search algorithm with n = 2; for more details see 
the main text. 



r P r (k = 1) P r (k = 2) Qr_ 

0.2500 0.5000 

1 0.5781 0.8750 1 

2 0.9249 0.9980 4 

3 0.9996 1.0000 13 

4 1.0000 1.0000 40 



For the case of k — 1 there is a single satisfying input 
\s). The probability of the algorithm succeeding depends 
on the projection of the final state onto the satisfying 
input and is given by 

P,H(s|K|00)| 2 (4) 

which simplifies to 

P r = 1 - (3/4) 3 ". (5) 

Clearly this converges rapidly to one, as shown by the 
numerical values listed in table [I] For the case k = 2 the 
target state is an equally weighted superposition of the 
two satisfying states, and the success probability is given 
by 

P r = 1 - (1/2) 3 " (6) 

which rises even more rapidly. 

Note, however, that r is the order of the recursive op- 
erator and not the number of queries, which is given by 
Q r = (3 r — l)/2. This query count, which is also listed 
in Table[IJ also rises rapidly with r, so the fixed-point al- 
gorithm is less efficient than a traditional Grover search 
when the value of k is known. This is unsurprising, as 
the traditional search is known to be optimal jj] in this 
case! When the value of k is unknown, however, the new 
algorithm outperforms previously known methods, as de- 
scribed by Grover 0. It is also relatively robust against 
experimental errors, as discussed below. 

III. EXPERIMENT 

This algorithm was implemented on a two qubit nu- 
clear magnetic resonance (NMR) quantum computer 
0- 0, ITol ITU IT^ . Il3| . searching for either one or two 
satisfying inputs for a function with four inputs. The 
spin system chosen was formed by the H and 13 C nu- 
clei in a sample of 10 mg of 13 C labeled sodium formate 
(Na + HCO^") dissolved in 0.75ml of D2O at a tempera- 
ture of 20°C. 

All experiments were performed on a Varian Unity 
Inova spectrometer with a nominal 1 H frequency of 
600 MHz. The 1 H and 13 C frequencies were adjusted to 



be in exact resonance with the respective nuclei so that 
the spin Hamiltonian in the rotating frame, written us- 
ing Product Operator notation Q , is given by the Ising 
coupling 

H = ttJ2H z C z (7) 

with J = 194.8 Hz. The measured relaxation times were 
T? = 6.5 s, Tf = 1.2 s, if = 16 s, and T^ = 0.6 s. A 
repetition delay of 120 s was used in all experiments; this 
more than seven times the longest Ti and so saturation 
effects can be ignored. The radio frequency (rf) pulse 
powers were adjusted so that a 90° rotation took 15 /is 
for both spins. 

The new search algorithm requires the implementation 
of U gates and R gates, as well as the inverse opera- 
tions. The U gates were implemented as simultaneous 
90° pulses, while the R gates were decomposed into peri- 
ods of evolution under the Ising coupling and composite 
^-rotations [l^, constructed from x and y-pulses. Rf 
gates were implemented for each of the four functions / 
with k = 1 and the six functions with k = 2; note that 
the Rq gate is identical to Rf for the function with 1 00) 
as the single satisfying input. Each gate was locally op- 
timized by combining pulses, but no optimization across 
gates was performed. To explore the effects of system- 
atic errors, experiments were performed u sing both naive 
rf pulses and BB1 composite pulses Il7| which cor- 
rect for systematic errors in pulse lengths arising from rf 
inhomogeneity. 

A pseudo-pure initial 1 00) state was prepared by spatial 
averaging |12j . Experiments were performed for each of 
the functions with k = 1 and k — 2, with the order r 
of the recursive search operator taking the values r = 
0, 1, 2, 3. The r — > 00 limit was simulated by directly 
transforming the initial state into the desired target state. 
The state of the spin system was then probed to obtain 
information on the performance of the algorithm. 

For the case of k = 1 the target state is a single eigen- 
state and the analysis is simple. A crush gradient was 
applied to dephase any off-diagonal error terms in the 
density matrix and the 1 H NMR spectrum was ob- 
served after a 90° 1 H pulse. The state of the first qubit, 
stored on the : H nucleus, is then encoded in the sign of 
the NMR resonance, and the state of the second qubit, 
stored on the 13 C nucleus, is revealed by which of the two 
components of the 1 H multiplet is observed. For more de- 
tails see ^3 ■ Note that this approach is only practical in 
systems with relatively small numbers of qubits, and in 
large spin systems it would be necessary to measure all 
the spin states directly. 

Finally, the success probability P r of the algorithm can 
be estimated from the intensity of the NMR spectrum 
compared with a reference spectrum . Because NMR 
experiments are only sensitive to the traceless part of the 
density matrix the observed signal strength is not directly 
proportional to P r , but rather to the fractional signal F r , 
which is given by F r = (4P r — l)/3. In particular, no 
signal is expected for the case P r — 1/4. 
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FIG. 1: Experimental 1 H NMR spectra from an implementation of Grover's fixed-point quantum search algorithm on a two 
qubit NMR quantum computer with one matching state. For a description of the readout scheme see the main text. Spectra 
are shown for the cases r = 0, 1, 2, 3, and a simulation of the r — > oo limit was obtained by directly transforming the initial 
state into the target state. Spectra are plotted using NMR conventions, with frequency (measured by the offset from the rf 
frequency) increasing from right to left. A horizontal axis is plotted below the bottom left spectrum and can be applied to all 
spectra. The vertical scale is arbitrary, but the same for all spectra. 



For the case of k = 2 the situation is slightly more 
complicated, as the target state is a superposition of the 
two matching states, but we chose to analyze the data in 
the same way. For four of the six functions this results in 
a 1 H NMR spectrum containing both components of the 
multiplet, with the result encoded in the signs of these 
two lines, while for the other two functions no signal is 
expected. For the four functions giving rise to visible 
signals the success probability and fractional signal are 
related by F r = 2P r — 1, and so no signal is expected for 
the case P r = 1/2. 

IV. RESULTS 

We began by implementing the four functions with 
k = 1 using naive rf pulses, with the results shown in 
Fig. ^ The spectra all have the form expected, show- 
ing one major component in each multiplet. A positive 
line indicates that the first qubit is in state |0) while a 
negative line indicates state |1). Signal in the left hand 
component indicates that the second qubit is in state |0) 
while the right hand component indicates state |1). The 
minor signals visible on the other component of each mul- 
tiplet, as well as the minor phase distortions visible in 
some spectra, can be ascribed to errors in the implemen- 
tation. 

As expected the signal intensity initially rises towards 
the limiting value, although it seems to fall slightly at 
r = 3. This point is explored in more detail in Fig. |3 
which compares the integrated intensity of the largest 
component with the theoretically expected values. Ini- 
tially the experimental data points lie quite close to the 



theoretical line, but for r — 3 the match is much less 
good. We originally ascribed this to the effects of errors 
in the pulse sequence, and in particular to the cumula- 
tive effects of pulse length errors, and sought to reduce 
these effects by using BB1 composite pulses The 
outline results of this approach are also shown in Fig. [3] 
(raw data not shown). 

While this approach did give slightly improved results 
for r < 2 (especially for the case r = where the search 
operator comprises a single 90° pulse) , it does not remove 
the drop in intensity seen at r = 3, which we now believe 
arises from errors in the implementation of the R gates. 
The errors are different for the four different functions 
(the error arising from noise in the experimental spectra 
was estimated by repetition and is much smaller than the 
scatter observed), consistent with this suggestion. 

The lack of improvement from the use of BB1 pulses 
may seem disappointing, but is in fact quite interesting 
in its own right. We have assumed that the U operator is 
implemented by a 90° pulse, but Grover's algorithm can 
be made to work with many different operators 0. In 
combination with the fact that Grover's new algorithm 
always moves towards the target state, this makes the al- 
gorithm intrinsically tolerant of pulse length errors Q ■ In 
fact the experimental spectra observed are of remarkably 
high quality, given that the case of r = 3 corresponds to 
the implementation of 26 two qubit gates (13 instances 
of Rf and 13 of Rq) and around 200 rf pulses. 

Finally we consider the situation when there are two 
matching states, that is k = 2. The search operators 
were implemented directly, rather than by applying two 
single-match operators in sequence, and are slightly sim- 
pler than for the case of a single target state (in some 




FIG. 2: Experimental success probability for the cases r — 0, 1, 2, 3 for the four possible functions with a single matching state. 
Results are shown for both simple single qubit gates implemented by naive rf pulses and for BB1 composite pulses. Fractional 
intensities are obtained by comparing the spectral intensities with those in the spectra obtained by directly transforming the 
initial state into the target state, and these are converted to probabilities as described in the main text. The theoretical result 
is shown as a smooth curve even though it is strictly only defined for integer values of r; experimental values are shown as 
squares, circles, diamonds and stars for target states of |00), 1 01) , 1 10) and |11) respectively. 



cases the function operators Rf do not require two qubit 
gates). There are six possible search operators Rf, all 
of which were implemented using both naive and BB1 
composite pulses, but here we concentrate on two cases: 
firstly where the target states are 1 00) and |01) (giving a 
H spectrum with both components of the multiplet pos- 
itive), and secondly where they are 1 10) and 1 01) (giving 
a spectrum with the lefthand component negative and 
the right hand component positive). 




The experimental results from these cases are summa- 
rized in Fig. where the intensity of each spectrum was 
obtained using either the sum or the difference of the in- 
tegrals of the two components as appropriate. We show 
results obtained using naive pulses, but as before the re- 
sults with BB1 pulses were very similar. As before the 
experimental data broadly follows the theoretical curve, 
and this time the results for r — 2 and r = 3 have almost 
the same intensity, as predicted by Table [I] This slight 
improvement may reflect the fact that the Rf gates are 
slightly simpler for k = 2 than for k = 1. 
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FIG. 3: Experimental success probability for the cases r = 
0, 1, 2, 3 for two of the six possible functions with two match- 
ing states. For further details see Fig. [5] Experimental values 
were obtained using naive pulses and are shown as circles for 
matching states of 1 00) and 1 01) and stars for matching states 
of 1 10) and |01). 



V. CONCLUSIONS 



In principle one could also study the cases of k = 3 and 
k = 4, but these are not particularly interesting. There 
is a correspondence between the functions with k and 
those with N — k satisfying inputs, and for quantum ora- 
cles (controlled phase gates) the operators differ only by 
global phases. Thus the case of k = 3 is indistinguishable 
from that of k = 1, as phase shifting three states in one 
direction is equivalent to shifting the fourth state in the 
opposite direction. The case of k — 4 is trivial, as apply- 
ing phase shifts to all the states is nothing more than a 
global phase, and thus the function operator corresponds 
to the identity operation. 



Our experimental results are broadly consistent with 
those predicted for an implementation of Grover's new 
fixed-point quantum search algorithm. The observed suc- 
cess probability initially rises with the order of the (re- 
cursively defined) search operator, although a fall-off is 
observed at r = 3, which we ascribe to experimental er- 
rors in the operators Rf and Rq. Despite these experi- 
mental errors, the results are remarkably good given the 
complexity of the pulse sequences involved. 

As predicted by Grover , the algorithm is remarkably 
robust to systematic errors which are equivalent for a 
gate and its inverse. This is largely true of the U gates, 
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which are implemented using rf pulses, and the use of 
BB1 gates to correct systematic errors has little effect 
except in the case r — 0. It is not true for the R gates, 
as our implementation of B) is somewhat different from 
that of R, and the errors in these two gates will not be 
equivalent. 

This error-tolerance property can in principle be used 
to develop methods for more general correction of sys- 



tematic errors jjj , but we do not address this point here. 

Acknowledgments 

We thank the UK EPSRC and BBSRC for financial 
support. 



[1] L. K. Grover, Phys. Rev. Lett. 79, 325 (1997). 

[2] L. K. Grover, Phys. Rev. Lett. 80, 4329 (1998). 

[3] C. H. Bennett and D. P. DiVincenzo, Nature (London) 

404, 247 (2000). 
[4] C. Zalka, Phys. Rev. A 60, 2746 (1999). 
[5] M. Boyer, G. Brassard, P. H0yer and A. Tapp, Fort. Der. 

Physik 46, 493 (1998). 
[6] J. A. Jones and M. Mosca, Phys. Rev. Lett. 83, 1050 

(1999). 

[7] L. K. Grover, A different kind of quantum search, quant- 
ph/0503205. 

[8] D. G. Cory, A. F. Fahmy and T. F. Havel, in Proceed- 
ings of PhysComp '96 (New England Complex Systems 
Institute, Cambridge MA, 1996). 

[9] D. G. Cory, A. F. Fahmy and T. F. Havel, Proc. Nat. 
Acad. Sci. USA 94, 1634 (1997). 
[10] N. A. Gershenfeld and I. L. Chuang, Science 275, 350 
(1997). 



[11] J. A. Jones and M. Mosca, J. Chem. Phys. 109, 1648. 
[12] J. A. Jones, Prog. NMR. Spectrosc. 38, 325 (2001). 
[13] L. M. K. Vanderspen and I. L. Chuang, Rev. Mod. Phys. 

76, 1037 (2004). 
[14] O. W. S0rensen, G. W. Eich, M. H. Levitt, G. Boden- 

hausen and R. R. Ernst, Prog. NMR. Spectrosc. 16, 163 

(1983). 

[15] R. Freeman, Spin Choreography (Spektrum, Oxford, 
1997). 

[16] S. Wimperis, J. Magn. Reson. A 109, 221 (1994). 

[17] H. K. Cummins, G. Llewellyn and J. A. Jones, Phys. 

Rev. A 67, 042308 (2003). 
[18] J. A. Jones, M. Mosca and R. H. Hansen, Nature 393, 

344 (1998). 

[19] M. S. Anwar, L. Xiao, A. J. Short, J. A. Jones, D. Blaz- 
ina, S. B. Duckett, and H. A. Carteret, Phys. Rev. A 71, 
032327 (2005). 



